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The IceCube detector, which is embedded in the glacial ice at the geographic South Pole, is the first neutrino 
telescope to comprise a volume of one cubic kilometer. The search for neutrinos of astrophysical origin is among 
the primary goals of IceCube. Point source candidates include Galactic objects such as supernova remnants 
(SNRs) as well as extragalactic objects such as Active Galactic Nuclei (AGN) and Gamma- Ray Bursts (GRBs). 
Offline and online searches for transient sources like GRBs and supernovae (SNe) are presented. Triggered 
searches use satellite measurements from Fermi, SWIFT and Konus. Complementary to the triggered offline 
search, an online neutrino multiplet selection allows IceCube to trigger a network of optical telescopes, which 
can then identify a possible electromagnetic counterpart. This allows to probe for mildly relativistic jets in SNe 
and hence to reveal the connection between GRBs, SNe and relativistic jets. Results from IceCube's triggered 
GRB search and a first limit on relativistic jets in SNe from the optical follow-up program are presented. 



1. Introduction 

Gamma-ray Bursts (GRBs) are prime candidates 
for the production of the highest energy cosmic rays 
because of the enormous ener gy that is released in 
such an event [Waxmanl fl995j (O(10 51 - 10 54 erg x 
f2/47r) in 7-rays, where f2 is the opening angle of a 
possible beamed emission). If the prime engine ac- 
celerates protons and electrons with similar efficien- 
cies this would be sufficient energy to account for 
the observed ultra high energy cosmic rays. The ob- 
served 7-rays would originate from high energy elec- 
tron synchrotron emission and inverse Compton scat- 
tering, while high energy neutrons would escape the 
fireball's magnetic field and later decay to protons, 
which would be responsible for the high energy cos- 
mic ray flux seen on Earth. The observation of keV 
to MeV 7-rays confirms the presence of high energy 
electrons in the fireball; however, because high en- 
ergy protons are deflected in inter-galactic and the 
Galactic magnetic fields no direct observation of pro- 
tons from GRBs is possible. Nevertheless, if high 
energy protons are present in the fireball along with 
high energy electrons it is reasonable to assume that 
pions will be produced through interactions of pro- 
tons with 7-rays from electron synchrotron radiation 
near the source, w hich would give rise to neutrinos. 
iGuetta et al.l [2004] give a detailed account of the ex- 
pected neutrino flux from such interactions. Here we 
present the most recent search for neutrinos in tem- 
poral and directional coincidence with reported GRBs 
done in IceCube. 

Recent observations imply a common physical origin 
of long GRBs and core-collapse s upernovae (CCSNe): 
a mas sive stellar explosion (see IWooslev and Bloom] 
20061 for a review) . According to the collapsar model 
MacFadven and Woo slcv 1999j . long GRBs (duration 



> 2 s) have their origin in the collapse of a massive, 
rapidly rotating star into a black hole surrounded by 
an accretion disk. Relativistic jets with Lorentz boost 
factors of 100-1000 form along the stellar axis. This 
GRB-SN connection gives rise to the idea that GRBs 



and SNe might have the jet signature in common and 
a certain fraction of core-collapse SNe host soft rela- 
tivistic jets. SN jets are suggested to be equally ener- 
getic and more baryon-rich, hence they are only mildly 
relativistic. Such soft relativistic jets would become 
stalled in the outer layers of the progenitor star, lead- 
ing to essentially full absorption of the electromag- 
netic radiation emitted by the jet and at the same 
time an efficient production of high-energy neu trinos 
[Razzaaue et aljfeOQSl lAndo and Beacomll2005j |. This 
motivates a search for neutrino emission from SNe, 
as ne utrinos would be able to escape from within the 
star. lAndo and Beacoml [2005j | present a detailed cal- 
culation of the expected neutrino flux. 
A dedicated search for a neutrino signal in coincidence 
with the observed X-ra y flash of SN 2008D has been 
conducted by IceCube Abbasi et al.ll20fTa | in order 
to test the soft jet scenario for CCSNe. Early SN de- 
tections, as in the case of SN 2008D, are very rare 
since X-ray telescopes have a limited field of view. 
However, neutrino telescopes cover half of the sky 
with high sensitivity at any time. If neutrinos pro- 
duced in soft relativistic SN jets are detected in real 
time, they can be used to trigg er follow-up observa- 
tions [Kowalski and Mohrll2007| . This is realized with 
an online analysis of neutrino data combined with an 
optical follow-up. This optical follow-up program is 
independent of satellite detections and thus comple- 
mentary to the triggered searches mentioned above. 
It is sensitive to transient objects, which are either 
7-dark or missed by 7-ray satellites. In addition to a 
gain in significance, the optical observations may al- 
low to identify the transient neutrino source, be it a 
SN, GRB or any other transient phenomenon produc- 
ing an optical signal. Hence it enables us to test the 
plausible hypothesis of a soft relativistic SN jet and 
sheds light on the connection between GRBs, SNe and 
relativistic jets. 

While up to now the triggered search is performed 
offline on an entire dataset (~ 1 year of data) with 
time consuming reconstructions on a large computer 
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cluster, the optical telescopes have to be triggered in 
real time and thus require an online analysis of the 
neutrino data. 



2. The IceCube Detector 



Th e IceCube neutrino telescope [Achterberg et all 
200G] has been under construction at the geographic 
South Pole since 2004 and was completed in the 
Antarctic summer of 2010/11. It is capable of detect- 
ing high energy neutrinos with energies of O(100)GeV 
and is most sensitive to muon neutrinos with energies 
in the TeV range and above. High-energy muon neu- 
trinos undergoing charged current interactions in the 
ice or the underlying rock produce muons in neutrino- 
nucleon interactions. The muon travels in a direction 
close to that of the neutrino and emits Cherenkov 
light. The deep ultra clear Antarctic ice is instru- 
mented with light sensors thus forming a Cherenkov 
particle detector. After its completion it comprises 
a volume of 1km 3 with 5160 digital optical mod- 
ules (DOMs) attached to 86 vertical strings at a 
depth of 1450 m to 2450 m. Each DOM consists of 
a 25 cm diameter Hamamatsu photomultiplier tube 
(PMT) and supporting hardware inside a glass pres- 
sure sphere. Here we present the offline analysis of 
data taken with the 40-string detector configuration 
(2008/04/05 to 2009/05/20) and the 59-string con- 
figuration (2009/05/20 to 2010/05/25) as well as the 
online analysis of the data taken from the start of the 
follow-up program on 2008/12/16 to 2009/12/31. In 
the following the deployment stages will be referred 
to as IC40 and IC59. 
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Figure 1: Neutrino spectra of the GRBs used in the 
triggered offline analysis. The thin lines represent the 
individual bursts while the solid thick line represents the 
sum of al l burs ts. Finally, the dashed line shows the 
IWaxmanl [200 3] prediction normalized to the number of 
GRBs observed. 



are taken preferentially from Fermi (GBM), Konus- 
Wind, Suzaku WAM, and Swift in this order. The 
7-ray spectra reported by the satellites were used to 
calculate the ne utrino spectra and fl ux as outlined in 



Appendix A of iGuetta et al.l |2004j . GCN does not 



always report values for all of the parameters used in 
the neutrino sp ectrum calculation. I n that case aver- 
age values from Abb asi et al.l [201 lbj | are used for the 
parameters not measured by the satellites. 



3.2. Offline Event Selection 



3. Triggered Offline Analysis 

The offline analysis is triggered by satellite detec- 
tions. It looks for neutrinos coincident in space and 
time with the observed 7-radiation. 



3.1. The GRB sample 

During the IC59 data taking period, 105 GRBs were 
observed in the northern sky and reported via the 
GRB Coordinates Network (GCN 1 ). Of those GRBs 
9 had to be removed, because IceCube was not tak- 
ing physics data. The GRB localization is taken from 
the satellite that has the smallest reported error. The 
start (Tg-kaj-t) and stop (T s |; p) times are taken by 
finding the earliest and latest time reported for 7- 
emission. The fluence, and 7-ray spectral parameters 



1 http : //gen . gsf c . nasa . gov 



The analysis presented here was designed to be 
sensitive to neutrino production from P7 interac- 
tions in the prompt phase of the fireball. To sep- 
arate sjj£njy__fnmi ^background a Boosted Decision 
Tree [Hoecker et al. 2007] was trained. The analysis 
was then optimized for discovery with respect to the 
Boosted Decision Tree score. The optimized value re- 
sulted in a final data sample of 85% atmospheric neu- 
trinos and 15% miss-reconstructed cosmic ray muons 
in the off time data sample (any events not within 
±2 hours of a GRB). An unbin ned maximum likeli- 
hood search [Abbasi et alJ[20Tfl | was performed and 
each event passing the boosted decision tree cut was 
assigned a probability of being a signal event from a 
GRB or a background event. 



3.3. Results from triggered offline 
analysis 

No events were found in the on-time data to be 
on-source (within 10° of a GRB) and on time with 
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Figure 2: The combined limit of the IC40+59 analysis is 
shown in add i tion t o th e flux prediction s from 
iGuetta et all |2004l ] and lWaxmanl [20031 ]. 



a GRB in the IC59 data set. In total 24 back- 
ground events (not necessarily on source) were ex- 
pected to be in the total time window and 21 were 
observ ed (none on-source). From the IGuetta et al.l 



20041 ] model 5.8 signal events were predicted and a 
final upper limit of 0.46 times the predicted flux can 
be set. This limit includes a 6% systematic uncer- 
tainty. The corresponding mo del dependent result 
presented in a previous analysis [Abbasi et al.ll2011b| 
sets a limit of 0.82 of the model flux. Combining both 
results in a limit of 0.22 t imes the flux calculated ac- 
cording to IGuetta et al.l [2004] . This result implies 
that the model in question is strongly disfavored. The 
caveat is that certain parameters in the model, such 
as the bulk Lorentz factor T, are not (or only rarely) 
measured and therefore theoretically calculated values 
with large uncertainties are used. A lower limit on the 
Lorent z factor is establish ed by pair production argu- 
ments [Guetta et al. I I2004JI. but the upper limit is less 
clear Recent papers Bregeon et al.ll201l[|Liang et all 
l2010LlSo"derberg and Ramirez-Ruizl l2003] suggest that 
T can take values of up to 1000 (316 was used in this 
analysis), but it is questionable whether this is valid 
for a large fraction of GRBs. Whether the uncertain- 
ties on other parameters can account for the strict 
upper limit is currently being investigated. 



4. Online Analysis 
Program 



Optical Follow-Up 



The online analysis triggers a network of optical 
telescopes by the detection of neutrino bursts (at least 
two neutrinos in a 100 s time window and an angu- 
lar difference between their reconstructed directions 
of less than 4°). 



4.1. Online Event Selection 

In order to rapidly trigger optical telescopes the first 
online analysis of high-energy neutrinos detected by 
IceCube was developed and implemented. The on- 
line analysis presented here was designed to be sen- 
sitive to neutrino production in soft relativistic SN 
jets. Straight cuts are applied in order to separate 
signal from background events. Restricting the search 
to the Northern hemisphere and imposing require- 
ments on the event reconstruction quality (e.g. the 
number of hits with small time residual or the like- 
lihood of the reconstruction) allows a suppression of 
the mis-reconstructed muon background. To suppress 
the background of atmospheric neutrinos, which we 
cannot distinguish from the soft SN neutrino spec- 
trum, we require the detection of at least two events 
within 100 s and an angular difference between their 
two reconstructed directions of A 1 ? < 4°. The choice 
of the time window size is motivated by the observed 
7-ray emission from long GRBs (typically 50 s) , which 
roughly corresponds to the time for a highly relativis- 
tic jet to penetrate the stellar envelope. The angular 
window A "I 1 is determined by the angular resolution 
of IceCube and was optimized along with the other se- 
lection parameters. The final set of selection cuts has 
been optimized in order to reach a multiplet rate of 
~ 25 per year corresponding to the maximal number 
of alerts accepted by ROTSE. The final data stream 
consists of 37% (70%) atmospheric neutrinos for IC40 
(IC59). Combining the neutrino measurement with 
the optical measurement allows the cuts to be relaxed 
yielding a larger background contamination and at the 
same time a higher signal passing rate. A doublet is 
not significant by itself, but may become significant 
when the optical information is added. The resolu- 
tion of the doublet direction is ~ 0.8°. 



4.2. Search for Optical Counterparts 

The IceCube multiplet alerts are forwarded to the 
robotic optical transient search experiment (ROTSE) , 
which consists of four identical telescopes each with 
a 45 cm mir row located in Aus tralia, Texas, Namibia 
and Turkey [Akerlof et al.| [2~003] . The telescopes stand 
out because of their large field of view of 1.85° x 1.85° 
and a rapid response. Once an IceCube alert is re- 
ceived by one of the telescopes, the corresponding re- 
gion of the night sky will be observed within seconds. 
Follow-up observations are performed regularly for 3 
weeks. In the initial phase with IC40 and IC59, the 
online processing latency of several hours made the 
search for an optical GRB afterglow unfeasible. We 
therefore focus on the SN light curve detection in the 
ROTSE data. 

Image subtraction followed by a detection algorithm 
to identify transient source candidates are applied to 
find an optical counterpart. 
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Figure 3: Limits on the c hoked jet SN model 
[Ando a nd Beacom 2 005l | for different boost Lorentz 
factors r as a function of the rate of SNe with jets p and 
the jet energy £j e t (colored regions are excluded at 90% 
CL). Horizontal dashed lines indicate a fraction of SNe 
with jets of 100%, 10% or 1%. 



boost Lorentz factors r = 6, 8, 10 as a function of the 
rate of SNe with jets p and the jet energy £j C t- Sys- 
tematic errors related to the simulated neutrino sensi- 
tivity and the SN sensitivity are included in the limit 
calculation. The 90% confidence regions for each re- 
value are displayed in the £j e t-/o-plane in Fig. |3] The 
most stringent limit can be set for high T-factors. A 
sub-population fo SN with typical values of T = 10 
and Ejct = 3 • 10 51 erg does not exceed 4.2% (at 90% 
confidence). This is the first limit on CCSN jets using 
neutrino information. 

Because of the successful operation of the optical 
follow-up program with ROTSE, the program was ex- 
tend ed in August 201 to the Palomar Transient Fac- 
tory Rau et ahl 12009} . which will provide deeper im- 
ages and a fast processing pipeline including a spec- 
troscopic follow-up of interesting SN candidates. Fur- 
thermore, an X-ray follow-up by the Swift satellite of 
the most significant multiplets has been set up and 
started operations in February 2011. 



4.3. Results of the Optical Follow-up 
Program 

Here the results from the analysis of data taking 
in the period of 2008/12/16 to 2009/12/31 are pre- 
sented. Table Q] shows the number of detected and 
expected doublets and triplets for the IC40 and the 
IC59 datasets as well as the number of detected and 
expected optical SN counterparts. The IceCube ex- 
pectation based on a background only hypothesis was 
obtained from scrambled datasets. The number of 
doublets shows a small excess, which corresponds to 
a 2.1 er effect and is thus not statistically significant. 
The expected number of randomly coincident SN de- 
tections, -/Vgjjj = 0.074, is based on an assumed core- 
collapse SN rate of 1 per year within a sphere of ra- 
dius lOMpc, i.e. 2.4 • 10^ 4 y^ 1 Mpc~ 3 , and a Gaus- 
sian absolute magnitude distribution with mean of 
— 18 mag and standard deviation of lmag for CCSN 
Richardson et al. 120061 ]. No optical SN counterpart 



was found in the data. 
We obtain the confidence level f or different combi- 
natio ns of SN model parameters [Ando and Beacom] 
l2005j] by using a pre-defined test statistic based on a 
likelihood function. The limit is calculated for the jet 



Table I Measured and Expected Multiplets 





SN 


Doublets 
IC40 IC59 


Triplets 
IC40 IC59 


measured 





15 19 





expected 


0.074 


8.55 15.66 


0.0028 0.0040 
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